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Abstract. Let p be a prime and a a quadratic non-residue mod p. Then the 
set of integral solutions of the diophantine equation x^ — ax\ — px\ + apx\ = 1 
form a cocompact discrete subgroup r Pitt C SL(2,R) and is commensurable 
with the group of units of an order in a quaternion algebra over Q. The problem 
addressed in this paper is an estimate for the traces of a set of generators for 
r Pl£l - Empirical results summarized in several tables show that the trace has 
\Q significant and irregular fluctuations which is reminiscent of the behavior of 

I the size of a generator for the solutions of Pell's equation. The geometry and 

arithmetic of the group of units of an order in a quaternion algebra play a key 
l i role in the development of the code for the purpose of this paper. 

H 

Introduction 

B The solutions of the Poll's option - *• _ 1 tau , discrete subgroup of au 

algebraic torus and is isomorphic to Z. An estimate of the form for the size of 

i— I a generator of this group is given by theorems of Schur and Siegel. This estimate 

reflects the "worst case" scenario and empirical results show that the actual size 
of the generator is often much smaller. Here we consider the diophantine equation 
Xq — ax\ — px\ + apx\ = 1, where a is a quadratic non-residue mod p. This 

^vq diophantine equation is analogous to Pell's equation in the sense that the integral 

solutions of this equation form an arithmetic subgroup r PiQ of SL(2,M) and is 
commensurable with the group of units of an order in a quaternion algebra over Q. 
But unlike the case of Pell's equation, this is a non-commutative (finitely generated) 
group and the question analogous to the size of the generator for Pell's equation is 
an estimate for xq (or equivalently the traces) for a set of generators for this group. 

• »— j In this note we obtain empirical estimates for these quantities for a set of p's and 

r^J a's, and according to the discriminant of the corresponding quaternion algebras. 

These results are displayed in Figures 1-4 in §5 which confirm the rather irregular 
behavior of the size and number of generators as a function of the discriminant. 

In order to find a set of generators for the group T p ^ a it is necessary to invoke the 
geometry of the fundamental domain for the action of the group on the Poincare 
disc. Once a fundamental domain is obtained a set of generators for the group can 
be exhibited and the main difficulty in the computer implementation is to make sure 
that a precise fundamental domain has been explicitly constructed. The key idea is 
due to Hutchinson (in a 1907 paper [Ht]) who introduced a very useful ordering of 
the group elements and a construction for a fundamental domain which was later 
popularized by Ford and is usually known after the latter. Hutchinson notes that 
his method is different from those in Frickc and Klein and lends itself much better 
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to numerical implementation and explicit constructions. Johansson noticed that 
by using only the geometry of a Ford fundamental domain (in the cocompact case) 
one can give an estimate for the size of generators for T p ^ a which is a cocompact 
subgroup in our case since a is a quadratic nonresidue mod p. 

The problem considered in this note has application to the explicit construc- 
tion of the fundamental group of a compact orientable surface as a group of 2 x 2 
matrices which is often encountered by students of mathematics. Poincare's geo- 
metric construction does not provide a discrete subgroup of PSL(2, R) sufficiently 
explicitly for displaying the desired matrix representation. This article in fact grew 
out of an attempt to provide a method amenable to a computer implementation 
for explicitly exhibiting sets of generators for the fundamental groups of surfaces 
defined as arithmetic quotients of the upper half plane H. 

In sections 1 and 2 basic relevant facts about quaternion algebras and their 
orders are recalled to fix notation and the geometry of the fundamental domain is 
discussed in section 3. An account of the relevant aspects of the works Hutchinson 
and Johansson is given in §4 and applied to the case under consideration in §5. 

It should be pointed that in a series of papers, J.H.H. Chalk [Ch-K & Ch-A] 
obtained estimates for the size of generators for solutions of quadratic forms based 
on methods of analytic number theory. However his estimates, when specialized to 
our case, are crude compared to those that one obtains by the method of Johansson. 
The method of Corrales et al [Co] has similarities to that of Johansson but the 
work of the latter gives better effective bounds. In [A-B] a method is described for 
explicitly exhibiting a set of generators for the group of units of Eichlcr orders in 
some special classes of quaternion algebras. 

1. Quaternion Algebras 

Definition 1. Let K be a field with characteristic ^ 2. H is a quaternion algebra 
over K if it is a four dimensional algebra over K with a basis {1, i, j, ij} satisfying 
the product rules 

i 2 = a, j 2 = b, ij = -ji, a,b £ K* . 

Then it is well known that H is a central simple i^-algebra, and is denoted 

by H = (jf^ ■ There is an involution h h taking h = a + ib + jc + ijd to 

h = a — ib — jc — ijd. The reduced trace and reduced norm of h £ H are defined as: 

tr(/i) = h + h, n(h) = hh 

By the Nocthcr-Skolcm Theorem for every separable quadratic if-subalgebra L 
of H we have H = L + Lui, for some lo £ H with following properties: 

u? =6 £ K* 

loI =Iuj, for all I £ L. 

We refer to this representation as H = {L, 9}. 

If E/K is any field extension, we have (^f) ®E= (^jf), and {L, 9} <g> E = 

{L ® E,6}. li H ® E = M(2, K) we say that E is a splitting field of H. One can 
prove that a quadratic field extension E/K splits a division quaternion algebra H 
if and only if E is iiT-isomorphic to a maximal subficld of H containing K. 

In this report, we consider the case where K is a number field with the ring of 
integers R. For every place v of K, let K v denote the completion of K at v. 
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Definition 2. Let H be a quaternion if -algebra. Then H v := K v <g> iJ is a 
quaternion i^- algebra. If i?„ is a division algebra, we say that if is ramified at v; 
otherwise we say H splits at v. 

Theorem 3. There are exactly two quaternion algebras over K v up to isomorphism; 
one is the matrix algebra M(2,K V ), and the other is M v = {L nr ,Tr}, where L nr is 
the unique unramified quadratic extension of K v , and ir is a uniformizer of K v . 

For H = (^^j and v a place of K, the Hassc invariant is defined as 

, _ J 1, if H splits at v ; 
£(a ' > v ~ \ -1, if if is ramified at v. 

If K = Q, the Hasse invariant coincides with the Hilbert symbol (a, b) v : 

, . _ J 1 , if X 2 — aY 2 — oZ 2 is isotropic over K v ; 
(a ' > v ~ \ -1, otherwise. 

For computing (a, the following properties of the Hilbert symbol are useful: 



(1) {a 1 bc) p = (a,b) p (a,c) p . 

(2) if p is a real place, then (a, b) p — — 1 iff a < 0, 6 < 0. In the case (a, 6) p = —1 
we say that H is definite, and otherwise H is indefinite. 

r 1, ifa,6ei?; 

(3) Ifp^ 2 and a and 6 are co-prime, then (a, = J (orp{a,p{6); 

[ (|), if p\a,p\b. 

(4) The number of places that is ramified in H is finite and Y\ p (a,b) p = 1, 
Va, b e K *, so that (a, 6) 2 = Ilp^2( a ' b )p 

Theorem 4. (1) Two quaternion K-algebras are isomorphic if and only if 
they are ramified at the same places. 
(2) Given an even number of non- complex places of K , there exist a quaternion 
K-algebra that ramifies exactly at these places. 

Proof. See [M-R] Theorem 7.3.6. □ 

Definition 5. The (reduced) discriminant dn of a quaternion algebra H/K is the 
product of prime ideals of R that ramify in H . We will refer to the set of ramified 
primes by Ram(iJ). 

By the Hilbert Reciprocity Law Ram(iJ) is finite of even cardinality. In the case 
of K = Q, dn is an integer, which if H is definite, is a product of an odd number 
of distinct primes, and if H is indefinite, is a product of an even number of distinct 
primes. In this paper, we mainly concerned with the indefinite cases. 

Let H = (^r) be an indefinite quaternion algebra, with a > 0. We fix the 
embedding $ : H — > M(2,Q v / a) with 

ci -n k( , ■, ■, ••-> ( ^o + ^i\/a x 2 + x 3 ^/a \ 

For an indefinite quaternion algebra H over Q denote the discriminant by 
dtf = d = piP2---P2m- By the Chinese Reminder Theorem, we can choose an 
integer a such that ( ) = — 1 for all pi > 2. Furthermore by Dirichlet's Theorem 
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on primes in an arithmetic progression one can assume a is a prime p such that 
p = 5 (mod 8) . For such a p, the Hilbert symbol: 

{p,d) Pi = (p,Pi) Pz = (j^J = -1, for all pi > 2. 

Nowp = 5( mod 8) implies (^j = —1. Since p = 1( mod 4) we have (/p~^) = ^' 
and therefore 



2m / \ 



Hence (jj^j ramifies at all primes dividing d, and is unramified at p, and at all other 

odd primes. So (^jr) has the same discriminant as H, and H = (i^^J — (^jj^)- 
We will make use of this representation of an indefinite rational quaternion algebras. 

2. Quaternion Orders 

Definition 6. An i?-order O in a quaternion algebra H/K is a complete i?-lattice 
in H which is also a ring with unity. 

By an integer of H we mean an element a € H such that tr(a) and n(a) are in 
R. Each element of an order is an integer in H. For an order O, we define the dual 
lattice as 

O* = {x G H; tr(xO) C i?} 
which is a complete i?-lattice in iJ. 

Definition 7. Let O be an order of H. The reduced discriminant of O, which 
denotes by d(0), is the ideal C R. 

One can prove that if O = R[xi, ...,x^], then n(0# ) 2 = Rdet[tr(xiXj)]ij, 
1 < i, j < 4. So for two orders O' C O, the discriminants satisfies 

d(O') - [O : C]d(C) 

and d(O) = d(0') implies that O = O' . We say that an order is maximal, if it is 
not properly contained in any other order. An Eichler order is the intersection of 
two maximal orders. 

Proposition 8. Let H = {L, 6} be a quaternion algebra over a local field K v , and 
p be the maximal ideal of R v . Let O be a maximal order in H. If H = M(2,K V ), 
then O is conjugate to M(2,R V ), and d(O) = R v . Otherwise H = L + Luj is a 
division algebra, O = Rl + Rl^ is unique, and d{0) = p. 

Let if be a number field. An i?-order O in H/K is a maximal order if and only 
if O p :— R p ® O is a maximal i? p -order for each place p of K. This follows from: 

o = h n e> 2 n o 3 n ... n o p n ... • 

From this, and the relation d{O p ) = d(0) p , we conclude that 

O is a maximal order d(0) = ] [ £ p = d(H). 

peRam(H) 

This gives a criterion for the characterization of maximal orders by their discrimi- 
nants. 
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3. Arithmetic Groups 

Let TC — {z 6 C ; ^s(z) > 0} be the Poincare half-plane with the hyperbolic 
metric 

, 2 dx 2 + dy 2 

ds = = . 

V 

Let O be an order in H = 0q) , and 

O 1 ={xe O; n{x) = 1} 

be the group of proper units of it. Under $ (see the image of O lies in 
PSL(2,R). One can show that this group is a discrete subgroup of PSL(2,M.) 
[M-R, Chapter 8]. PSL(2,R) acts on 7i by Mobius transformations 

*■«-(: i) 6Psift,) - 

Then T := also acts on Ti properly discontinuous and freely i.e. has no 

fixed points. Since we can transform TL to the unit disc IA by the action of matrix 

ijj = I j J , which sends the i to the origin, the conjugation of T by ip gives an 
embedding f of O 1 in 5C/(1, 1). 

Definition 9. Given 7 = ( ° ^ e Si(2,C) such that c ^ 0, the circle C 7 := 

{zeC; |cz + rf| = 1} is called the isometric circle of 7. We denote by r 1 and o 1 
its radius and center, respectively. 

The isometric circle is characterized by the property that 7 acts on it as a 
Euclidean isometry. By a straightforward computation, the center and the radius 
of C 7 are the numbers o 7 = — d/c and r 7 = l/|c|, respectively. 

The action of PSL(2,M.) extends to TC U K U {00}, and we have three types of 
elements in the group action, depending on their fixed points. Fixed points are 
determined by the solutions of 

f(z) =cz 2 + (d-a)z-b = 0. 

We have three cases according to the sign of the discriminant 8(f) = (a + d) 2 — 4. 

Definition 10. Assume 7=1? ^ I defines a transformation different from 

\b c ) 

±Id. 7 is 

• hyperbolic, if it has two distinct fixed points in R U 00 , equivalently if 
|tr( 7 )| > 2; 

• elliptic, if it has two complex conjugate fixed points , equivalently if 

|tr( 7 )l <2; ' 

• parabolic, if it has a unique fixed point in I U 00 , equivalently if 
|tr( 7 )| = 2. 

By using the Hasse-Minkowski Principle on quadratic forms, one can prove that 
the quaternion groups, i.e. groups arising from the proper units of an order in an 
indefinite quaternion algebra have no parabolic elements. Elliptic and hyperbolic 

elements are diagonalizable as f ^ A^ 1 )' em P^ c elements, A 2 = e l6 , and 
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they are torsion elements, and for hyperbolic elements 7 7^ 1 is real. In the case of 
quaternion groups elliptic elements are only of orders 4 or 6, 

Definition 11. A connected closed polygon PCWURU {00} is a fundamental 
domain for the action of T on 7i if any two points z,z' in the interior of T> are not 
T-equivalent, i.e. there is not any 7 <G T such that z = 7(2'), and if each point in 
H is T-equivalent to some point in V. 

Ford describes a method to construct a fundamental domain for the action of 
any discrete subgroup of PSL{2, M) on the unit disc by using their isometric circles 
[Fl]. His method is applicable to quaternion groups of interest in this paper and 
the resulting fundamental domains have finitely many sides. In Ford's method 
the region which is exterior to all isometric circles of elements of T constitute a 
fundamental domain of the action of T on the unit disc. The set of elements 7 
corresponding to the sides of this region give us a set of generators of T; however 
this set is not necessarily minimal. 

By attaching a Dedckind zeta function to a maximal order O in an indefinite 
quaternion algebra, and relating the residue of its pole at s = 1 to the hyperbolic 
area of the fundamental domain of the action of the corresponding unite group T 
on 7i, Eichler [E2] showed that 

(3.1) Voi(r\w) = | H(p-i) 

p\d(0) 

Vigncras [V] recovered this result by using local-global Tamagawa measures and 
adclic techniques, and computed the area of the fundamental domain for a maxi- 
mal order with respect to the Poincare arithmetic metric ds := ds/2n. With this 
normalization the volume is the rational number | Yl p \d(o)(P~ -0 = ^(diO)). (In 
[E2] and [V] these computations are considered for the wider class of Eichler orders 
of level N) 

4. Fundamental Domains And Generators 

For finding the Ford fundamental domain and a set of generators, one computes 
isometric circles for elements T. For a computer implementation of this, it is nec- 
essary to order elements of f . Then isometric circles are computed one at a time 
according to this order until a closed domain 2) is obtained. This 2) is not neces- 
sarily a fundamental domain for T\U, and the actual fundamental domain may be 
smaller since the isometric circles of the other elements of f may cut across £>. To 
address these problems we proceed as follows. 

An upper bound for the generators is obtained in [Ch2] and [Ch-K] . Defining the 
norm of 7 e T C SLS(2, R) as 

|| T ||^a 2 + fo 2 + c 2 + d 2 for 7 =(" J), 

and assuming that T\H is compact, Chalk proved that there exists a set of gener- 
ators {A-l, ...,An} for T with the following properties: 

(1) PiH < ||A 2 || < ... < \\A N \\ 

(2) H^iH <ciJV 

(3) H^+iH < C2 ^VP*|| 5 (i = l JV-1). 
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He also computes the explicit constants c\ and ci in the case of maximal orders of an 
indefinite quaternion algebra over Q [Ch-K] , and obtained the the following bounds 
for the generators of the group of proper units of maximal orders in a quaternion 



algebra with the representation H 

N<Q + (t>{d H ), 
Pi|| < 2 + 4 7 r- 1 iV, 

The bound for \\Ai\\ is improved in [Ch2] and [Ch-A]. The exponent 5 for ||A; + i|| 
in terms of \\Ai\\ yields a crude and practically inapplicable estimate. 

For a faster implementation of the construction of the true fundamental region 
T> we use two methods described by Hutchinson [Ht] and Corrales et.al.[Co]. Let 
F(z) = 1 — zz, then F > represents the interior of the unit disc. One observes 
that if \z\ < \z'\, then l-FX^)! > |F(z')|. Therefore among the points of isometric 

circle C 7 associated to any 7 = ^ " 5 ) ^ ^' ^ assumcs a max ™ m on the 

nearest point to the origin which we denote by z 1 . The boundary of V consist of 
arcs of isometric circles that are closest to the origin and F assumes its maxima on 
these circles. We have: 

So the boundary circles for V occur among those 7's determined by the smallest 
values of |o|. Hence, we start with the small values for |a|, and compute the required 
norm for c by \c\ = \/\a\ 2 — 1, and the above descriptions assure us that the first 
time we attain a closed region, it will be a fundamental domain for f . 

The method of Corrales et.al. is based on the assumption that T is cocompact. 
If c G H is not a fixed point of T, and p denotes the hyperbolic distance, define 

D g (c) :={xeH; p(x, c) < p(x, g(c))}, geT 

to be the half plane containing c. Then 

Dr(c):= p| D g (c) 

is known as a Dirichlet fundamental domain of T\H. Now for a positive real number 
r set 

(4.2) D r (c) = f]{D g (c) I 1 + g e T with g(c) G B(c, r)} 

where B(c,r) = {x G TL ; p(x, c) < r}. Since -Dr(c) is compact, there exist a 
positive number R such that -Dr(c) = Dr{c). For finding such R we let k\, k 2l k 3 , ... 
be an increasing sequence of positive numbers (goes to infinity) , and compute 's 
until Dk n is compact for some fc„. By the relation [B, Theorem 4.2.1]: 

(4.3) ||0|| = 2coBhp(i, ff (t)) 

the computation of D^'s is reduced to those elements g G T for which \\g\\ < ki. 
Now set R = 2(max({fc„/2} U {p(c, x) ; x G D kn (c)})). For this R we see that 
^ D kn {c) C B(c,R/2). Also for any g G T which p(c,g(c)) > R, we have 
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B(c,R/2) C D g (c), because if there is some x £ B(c, R/2)\D g (c), then p(x,g(c)) < 
p( x ,c) < R/2, and so R < p(c,g(c)) < p(c,x) + p(x,g(c) < R, a contradiction. It 
follows that for the computation of the fundamental domain and the determination 
of generators, it suffices to consider elements g £ T with entries of g bounded by R. 

The method [Co] is computationally intensive since it requires the calculation 
of all the domains D^. . By a good choice of k\ the computational burden can be 
reduced. This method can be translate to the Ford fundamental domain, if instead 
of (|4.3|) we use the relation: 



(4-4) = ] , for 7 £ f . 

VIItII - 2 

Since T> is compact, then there is a disc 

U s = {z£C ; \z\ < 1-e}, <£ < 1, 

such that 2? C or equivalently for a set of generators {74} we have r~ li > e. So 
any good e which estimates r^. 's, gives an upper bound for entries of 7i's. Since we 
want that T> C U e , a crude bound can be assumed with respect to the area of T> as: 



k \ y 2 + Vol(X>) / 

where fc is larger than 2. In order to apply this method, it is necessary to know 
the area of T>. In the next section we introduce a method for computing Vol(2?) for 
quaternion algebras over Q. 

5. The Main Illustration 

Here we want to work in a special class of examples, and apply the above methods 
to obtain a bound on the trace and norm of a set of generators. We consider integer 
solutions of the quadratic Diophantine equation 

2 2 2 i 2 i 

x — ax 1 — px 2 + ap%3 = 1 

where p is an odd prime, and a < p is a quadratic non-residue mod p. Let 

x + xiy/a t/p(x2 + x 3y /a) \ i Xq — ax\ — px\ + apx\ = 1 
y/p(x 2 - xz^Ja) Xq~Xi^/cl ) x i £ Z Vi = 0, 3 

Then r p a is the group of units of the canonical order O p . a = Z[l,i,j, ij] in the 
quaternion algebra H — (^fj , and Xq — axf — px\ + apx| is the norm on H. 

Lemma 12. For p = 1 (Viod 4^ i/ie group T Pta is torsion free. 

Proof. Let Id^ 7 € r p Q be a torsion clement. Then its eigenvalues are e ±l0 and 
the characteristic equation is A 2 — (2cos#)A +1 = 0. Therefore xq — 0,±1, and 
7 is conjugate to a rotation matrix. In the latter case, xq — ±1, 7 is conjugate 
to ±Id, and so 7 = ±Id. In the former case, Xq = 0, the norm form becomes 
1 + ax 1 + px 2 = apx\ . Reducing mod p gives us 

1 + ax\ = mod p 

which is not possible, since —1 is a quadratic residue and a is not. □ 



r — 

L p,a — 
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Remark 13. The trace of an element of T pa is an even integer. Since an element 
of finite order is conjugate to a rotation matrix, every non-trivial element of finite 
order in r p a necessarily has trace zero. Let T p a be the subgroup of r pa , consisting 
of elements with xi = mod a. It is trivial that such matrices form a subgroup 
of finite index in r pa . For 7 S T® a we have £0 = 1 mod a, and therefore r p is 
torsion free. 

We apply Hutchinson's method to the given r PjQ . First we transform r Pja to the 
group Tp tU which acts on the unit disc: 

_ / x + ix 3y /ap xxsfa-ixzJp \ GTTM n 
7 = x + 1x1 + jx 2 + i]x 3 h-> 7 = ^ v . v 6 5(7(1,1) 

Since the Hutchinson's method works step by step to find the fundamental domain, 
we wrote a computer program which generates 7 G r p a , and stops when a compact 
domain was obtained. The following four diagrams are the bounds founded by the 
program experimentally for xqS, and number of generators, for some pairs (p,a). 
It is clear that the bounds for x\,X2 and X3 are comparable to those for xq. 



For these quaternion algebras, we can also apply the method of Johansson to find 
good bounds for Xi's. In order to do this we compute the area of the fundamental 
domain of the group of proper units A4 1 of a maximal or Eichler order M. D O p ^ a 
on U, and then apply the formula [V] : 

Vol(r Pl0 \«) = [M 1 : r Pi0 ] VolCM^W) 

to determine the area of the r p , a - Then (4.4 1 and (4.5) give bounds for the entries 
of 7's. The main problem here, is to compute the index [A^ 1 : r p a ]. There are 
two distinct methods for this purpose. One is given in [Chi]. This method is based 
on the reducing problem to the theorems about quadratic forms on finite fields by 
using [Chi, Theorem 2] 

[M 1 : r p .J = [] [(^V = FpM 

q prime, q\2m 

where m is an integer such that mM C O p ^ a , n = v q (m) is the g~adic valuation of 
m if q is odd, and n — max(8, v^irn)) when q = 2, and 

(r p , )«« = {u e Ptaq ; n{u) = 1( mod q n )}, 

(M 1 )^ = {veM q ; n(v) = 1( mod q n )}. 

For the computation of index, we combine the method described in [J3], the 
relation [V, IV.1.7]: 

[O 1 : O' 1 } = \\_pl : a''} which & C O, 

p 

where O p = O ® R p , and the Eichler invariant in the localizations of the orders. 
The Eichler invariant for O p ^ M(2,R p ) is defined as 

C 1 if O p /J(O p ) ^k p O k p ; 
e(O p ) = I if O p /J(O p ) S k p ; 

y —1 if O p /J(O p ) is a quadratic extension of k p , 

where k p is the residue field of K p and J(O p ) is the Jacobson radical of O p . These 
are the only possibilities for O p / J(O p ), because the dimension of semisimple algebra 
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Figure 1 . Number of generators plotted against the pairs (p, a) 




FIGURE 2. Bounds for xq's plotted against the pairs (p, a) 



GENERATORS OF ARITHMETIC QUATERNION GROUPS AND. 



11 





Figure 4. Bounds for Xq's ordered by discriminant dn 
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O p /J(0 p) over k p is at most 4. In fact, it cannot be a cubic or quartic field over 
k p because this would be a simple separable extension of k p , and it would follow 
from Hcnscls lemma that this can be lifted to a field contained in the algebra over 
the K p which is impossible. Also it cannot be a quadratic field times k p because 
idcmpotents can also be lifted, and then the algebra over K p would contain a 
corresponding subalgebra which is again impossible by the Double Commutator 
Theorem. So the only possibilities are those listed above. 

Lemma 14. In the rational quaternion algebra H p = H ® Q p , let two orders 
Op C M p be such that M p is maximal in H p and d(O p ) = p n . Then 

[Ml : ol] = [m; : o;]/[z; ■. n {o;)\ 

and 

p n-i {p 2 _ i)(p_ e(Op))- 1 if H p = M(2, v „ ; . 

l p . 



lM * . m _ J P n - 1 (p 2 -l){v-e{Op))- 1 if Hp 

l^p -U p l-^ p n-l {p+1){p _ e{0p)) -l ifHp 



Proof. See [Kr] Theorem 1. □ 

For an order O with the dual lattice C # and the dual basis {1,61,62,63} we 
assign the quadratic form 

fo{X 1 ,X 2 ,X 3 ) := d(0) • n(Xi6i + X 2 b 2 + X 3 b 3 ). 

Proposition 15. Let O be an order in the rational quaternion algebra H , and p 
be an odd prime. If fo = X\ + up r X 2 + wp s X^ with r < s and (uw,p) = I, then 

e(O p ) = l ^ r = 0, s > 1, (^-^j = !> 

e(O p ) = -l ^ r = 0, s > 1, (^^j = -1. 
In this case e(0 2 ) = 0. 

Proof. Sec Theorem 4.3 and Propositions 5.8 & 5.9 in [J3], and Satz 10 in [El]. □ 

Finally Johansson obtained the following formula for the volume of a fundamen- 
tal domain in the case of rational indefinite quaternion algebras [Jl, Prop. 5. 10]: 

q\d H q\d(0) 



■\d{0) J] [Z; : N{0* q ) 



q 2 



q ( q -e(O q) ) 

which M is the maximal order of H and q is prime. In the case of O x p a — T pM we 
have: 

fo V:a = -Xl + aX 2 2 +pXl 

Also for q an odd prime [Z* : N(0*)] = 1 since p is a prime and a < p (cf. 
[Jl] Lemma 5.3). For q — 2 the index depends only on the power of 2 in the 
decomposition of a. In fact, [Z 2 : N{0 2 )\ = 2 or 1 according as 4|a or 4 { a, 
respectively. 

In order to evaluate the efficacy of Johansson's method, we compared the results 
with the exact answers one obtains from Hutchinson's method. The results are 
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Figure 5. A comparison between the results of Hutchinson's 
method (the blue bars) and Johansson's bounds (the red bars), 
in logarithmic scale, about generators norm. 
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shown in Figure 5 (it is plotted in logarithmic scale). Note that Johansson's method 
can over-estimate or under-estimate by factors as large as 55 and 27 respectively, 
in the range of computed values. The average discrepancy with the exact values is 
almost a factor of 4. 

6. Open Problems 
In this section we briefly describe three related open problems. 

(1) One can similarly consider families of cocompact arithmetically defined 
discrete subgroup of SL(2, C) acting properly discontinuously on the upper 
half space and ask analogous questions about the size of a set of generators 
for such groups. 

(2) The results in this note have been empirical and we do not know of an 
analogue of theorems of Schur and Siegel for units of quaternion algebras 
which are definite over Q. It is possible that analytical methods based on L- 
functions and representation theory can be used to obtain better estimates 
for the size of a set of generators. 

(3) Here we only considered discrete groups commensurable with the group of 
units of a quaternion algebra definite over Q. One may ask similar questions 
about subgroups of finite index in groups of units of orders in quaternion 
algebras definite over a totally real number field and acting properly dis- 
continuously on a product of upper half planes. It is an interesting problem 
whether the construction of [A-B] for certain Eichler orders can be extended 
to this case. 

Acknowledgement - The author would like to thank Professors P. Bayer, Stefan 
Lemurell and Kleinert for very helpful communications, and especially Professor 
M. Shahshahani for suggesting the problem and consistent encouragement. 
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